1 



00 
On 

Or 
C/3 



Universal and non-universal behavior in Dirac spectra 
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We have computed ensembles of complete spectra of the staggered Dirac operator using four-dimensional SU(2) 
gauge fields, both in the quenched approximation and with dynamical fermions. To identify universal features 
in the Dirac spectrum, we compare the lattice data with predictions from chiral random matrix theory for the 
distribution of the low-lying eigenvalues. Good agreement is found up to some limiting energy, the so-called 
Thouless energy, above which random matrix theory no longer applies. We determine the dependence of the 
Thouless energy on the simulation parameters using the scalar susceptibility and the number variance. 
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1. Introduction 

The low-energy particle spectrum implies that 
the chiral symmetry which the QCD Lagrangian 
possesses in the limit of massless quarks is spon- 
taneously broken. The corresponding order pa- 
rameter, the chiral condensate, can be related to 
the pion decay constant through the Gell-Mann- 
Oakes-Renner relation [Q. On the other hand, 
the chiral condensate is directly related to the 
density of the smallest eigenvalues of the Dirac 
operator via the Banks-Casher relation 0. 

On the lattice, the density of the low-lying 
eigenvalues of the Dirac operator is strongly de- 
pendent on the choice of the fermion action. The 
staggered fermion action respects a chiral Ua(1) 
symmetry and, as a consequence, the nonzero 
eigenvalues of the Dirac operator occur in pairs 
±A„. Therefore, the eigenvalues near zero and 
in the bulk of the spectrum can be distinguished. 
The Wilson fermion action breaks chiral symme- 
try explicitly, and there is no Ua(1) symmetry for 
the hermitian Wilson Dirac operator. 

In this contribution, we will be concerned with 
universal features in the spectrum of the stag- 
gered Dirac operator. We will mainly concentrate 
on the low-lying eigenvalues since they are partic- 
ularly relevant for chiral symmetry breaking. The 
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Banks-Casher relation 7rp(0) = US, relates 
the eigenvalue density, p(X) — (J2 n ^(^ — ^n))> of 
the Dirac operator at zero virtuality to the abso- 
lute value of the chiral condensate, S. V is the 
four-volume. If chiral symmetry is spontaneously 
broken, this relation implies that the spacing of 
the small eigenvalues is ~ 1/(US). Since the 
eigenvalues of the non-interacting Dirac operator 
are spaced like I/U 1 / 4 , this means that the Dirac 
eigenvalues in QCD must be strongly correlated. 

Leutwyler and Smilga || derived an effective 
low-energy theory which is valid in the range 
1/A < L < 1/m-n, where L is the linear extent 
of the box, A is a typical hadronic scale, and rn^ 
is the pion mass. In this region, the kinetic terms 
in the chiral Lagrangian can be neglected, and 
only the symmetries are important. It was then 
conjectured by Shuryak and Verbaarschot j| that 
the distribution of the small Dirac eigenvalues, 
averaged over gauge field configurations, is uni- 
versal in the sense that it only depends on global 
symmetries of the Dirac operator. The essential 
ingredient is the spontaneous breaking of chiral 
symmetry in the QCD vacuum. Given this fact, 
the statistical properties of the low- lying eigenval- 
ues can be computed in a much simpler theory, 
e.g., chiral random matrix theory (RMT). Alter- 
natively, one can use the finite-volume partition 
function of Leutwyler and Smilga. We will not 
address this point since it is discussed in the con- 



2 



tribution by P.H. Damgaard g. 

To resolve the low-lying Dirac eigenvalues, one 
rescales the energies by a factor of V£ and defines 
the microscopic spectral density m, 



topological charge v is given by 



(1) 



^ (z) = y lhn oo^(^) • 

This is a typical universal quantity which can be 
computed in chiral RMT. Other universal quanti- 
ties are the distribution of the smallest eigenvalue 
and higher order spectral correlation functions on 
the microscopic scale A ~ 1/(VS). 

We note that by concentrating on the low-lying 
eigenvalues, we are considering a finite- volume ef- 
fect. In particular, by computing analytical re- 
sults we hope to learn something about the ap- 
proach to the thermodynamic limit. On the lat- 
tice, one necessarily has to perform an extrapola- 
tion to this limit. An example of the utility of the 
RMT results with regard to the thermodynamic 
limit has already been given in Ref. . 

We briefly review chiral random matrix theory 
in Sec. |^ and describe numerical details of our 
simulations in Sec. ||. Lattice data for the small 
Dirac eigenvalues obtained in the quenched ap- 
proximation are compared with RMT predictions 
in Sec. ^. New results with dynamical fermions 
are discussed in Sec. [|. These are particularly 
interesting with regard to the chiral limit. Of 
course, the random-matrix approach only works 
in a limited energy range. In Sec. ^, we will iden- 
tify the domain of validity quantitatively. We 
conclude with a summary and an outlook to fu- 
ture work in Sec. 

2. Chiral random matrix theory 

The Dirac operator in the continuum is defined 
as D = 7 M (9 M +gAfj), where A denotes the gauge 
field and g is the coupling constant. This opera- 
tor is anti-hermitian. Because of the Ua(1) sym- 
metry {75, D} = 0, all nonzero eigenvalues of iD 
come in pairs ±A„. In a chiral basis, the Dirac 
matrix has the structure 



iD = 



T 
Tt 



(2) 



where T is some complicated matrix. The Eu- 
clidean QCD partition function in a sector of 



Z, 



QCD 



= / VA^e- s & Y[ det ( D + 771 f) 
J f=i 

VA^e-^Y[m^l[(Xl+ m }), (3) 



A„>0 



where S s \ is the gluonic action, Nf is the num- 
ber of quark flavors with masses m/, and the 
A„ are the eigenvalues of iD. The superscript 
(y) on A means that the path integral is only 
over gauge fields with fixed topological charge v. 
The partition function is then given by -Zqcd = 

J2„ e^-^QCD' wnere is the vacuum angle. 

The basic idea of the random matrix approach 
is to substitute the matrix T in Eq. (|^) by a ran- 
dom matrix W, respecting the global symmetries 
of the problem. If the dimension of the matrix 
W is N x (JV + v), there are v exact zero modes. 
The RMT partition function reads 



7 (y) 

J RMT 



JV, 



VWP {W) Y[ AetiWW^ +m 2 f ), (4) 



/=i 



where the precise form of the distribution, Po(W), 
replacing e _5gl is unimportant, provided that it 
is invariant under rotations of W and that we are 
interested in universal properties |Q. For con- 
venience, one often uses a Gaussian distribution. 
Eventually, one is interested in the limit N — > 00 
which can be identified with the thermodynamic 
limit. The classification of the random matrix en- 
sembles pertaining to different symmetry classes 
can be found in Ref. || . We will discuss staggered 
fermions with gauge group SU(2) which are de- 
scribed by the chiral symplectic ensemble (chSE) 
of RMT. In this ensemble, the elements of W are 
real quaternions. Diagonalizing W and express- 
ing it in terms of angles and radial coordinates, 
the distribution of W can be written as || 

P(W) = P (W) Y[ deiiWW^ + m)) 



P Q ({\}) A 4 (A 2 ) J] A^ l+3 II ( A « + ™/) . (5) 

n /=! 
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where A(x) = n^j^i — x j) 1S the Vandermonde 
determinant. The main mathematical problem 
then consists in performing integrations of this 
expression over all but a few variables A„ in the 
limit N — ► co. 

3. Numerical details 

To check RMT predictions, one often replaces 
an ensemble average by a spectral average which 
is possible due to spectral ergodicity. Here, how- 
ever, we are interested in the low-lying eigenval- 
ues, and a spectral average is not possible. Thus, 
we need good statistics. 

Recently, it has become feasible to calculate 
large ensembles of complete spectra of the stag- 
gered Dirac operator in SU(2) Q on lattices of 
size up to 16 4 , extending the methods developed 
by Kalkreuter 0. We deter mine the complete 
spectrum in order to make sure that the distri- 
bution of the low-lying eigenvalues has no other 
numerical uncertainties than the finite precision 
in 64 bit arithmetic. 

Our numerical simulations were done on a 
CRAY T3E. The configurations were generated 
using a hybrid Monte Carlo algorithm. The 
boundary conditions were periodic for the gauge 
fields and periodic in space and anti-periodic in 
Euclidean time for the fermions. For the diago- 
nalization of the Dirac matrix, we employed the 
Cullum-Willoughby version of the Lanczos algo- 
rithm [fTo|| . In SU(2) with staggered fermions, ev- 
ery eigenvalue of iD is twofold degenerate because 
of a global charge conjugation symmetry. In ad- 
dition, the squared Dirac operator —D 2 couples 
only even to even and odd to odd lattice sites, re- 
spectively. Thus, on a lattice with N sites, — D 2 
has N/2 distinct eigenvalues. There is an exact 
sum rule for the distinct eigenvalues of —D 2 , 

(Ka) 2 =N, (6) 

A n >0 

where a is the lattice constant. Since we gener- 
ated complete spectra, we could use this sum rule 
to check the accuracy of our eigenvalues. Equa- 
tion (^) was satisfied with a relative precision 
of about 1CT 8 . In Table |, we summarize our 
quenched spectra. As an example, for L = 16 the 
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0.00401(3) 
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0.0293(3) 
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2.2 
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0.0089(1) 


1.2(2) 


2.4 


16 
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0.00390(9) 


1.2(3) 


2.5 
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576 


0.194(9) 


8(3) 


2.5 


16 


543 


0.016(2) 


12(7) 


Table 1 



Summary of our quenched spectra with j3 — 4/ g 
and V = L 4 . The last two columns represent the 
average value of the smallest eigenvalue (in units 
of l/2a) and its integrated autocorrelation time, 
respectively. 



eigenvalue spectrum has 32,768 different elements 
spanning four orders of magnitude. From Fig. |^ 
below one can read off that the fluctuations of the 
smallest eigenvalue cover one order of magnitude. 



4. Quenched results 

Our first studies were done in the quenched ap- 
proximation where the fermion determinants in 
Eq. (|B|) are absent. Most of the results were re- 
ported in Refs. §,§,[ll|- In Fig. |, we show the 
microscopic spectral density obtained for three 
different sets of parameters, covering strong and 
weak coupling. To compare the data with the 
RMT prediction [Eq. (@) below with a = 0] 
one needs to compute the parameter VY, which 
sets the energy scale. This parameter is deter- 
mined by the data via the Banks-Casher rela- 
tion, VTi — 7rp(0). [Since we have the complete 
spectrum for many configurations, we can eas- 
ily fit jo(0)]. Thus, the comparison of Fig. |l| is 
parameter-free. 

The microscopic spectral density has an oscilla- 
tory structure. The maxima correspond roughly 
to the most likely positions of the individual 
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Figure 1. Microscopic spectral density of the lattice Dirac operator (histograms) and RMT prediction 
(dashed lines) for three different parameter sets. Lattice size L 4 , inverse coupling constant j3 = 4/g 2 , and 
number of configurations are shown in the figures. 



eigenvalues. One can also compute the distribu- 
tion of the smallest eigenvalue alone, or higher 
order spectral correlation functions. For details, 
we refer to Refs. §,0. 

The data agree with the RMT prediction for 
topological charge v = (an example is shown 
in Fig. |f). This is due to the fact that in the 
derivation of the RMT results for v ^ 0, it is 
assumed that the Dirac operator has \v\ exact 




0.01 



Figure 2. Distribution of the smallest eigenvalue 
of the lattice Dirac operator (histogram) and 
RMT predictions for topological sectors v = 
(short dashes) and v = 1 (long dashes) for (3 — 2.4 
on a 16 4 lattice. 



zero modes. On the lattice, the would-be zero 
modes (in lattice units) are shifted by an amount 
proportional to a 2 . Therefore, unless a is very 
small or some form of improvement is used, we 
expect to find agreement of the lattice data with 
the RMT results for v = 0. 

We can see from Fig. |l| that the agreement of 
the data with the RMT result breaks down at 
some value of the rescaled energy. The domain 
of validity of RMT depends on the lattice volume 
and on the coupling constant. A quantitative dis- 
cussion of this issue will be given in Sec. ^. From 
the figures, one can see that the range of valid- 
ity increases with increasing lattice volume and 
decreasing [3. 

5. Results with dynamical fermions 

The results obtained in the quenched approx- 
imation are very encouraging. The natural next 
step is to include the fermion determinants in 
Eq. (H|). Now an additional parameter enters 
the problem, the quark mass (assuming that we 
consider degenerate quarks). Again, we are in- 
terested in the small eigenvalues of magnitude 
~ 1/(7E), i.e., z = XVY, ~ 0(1). From Eq. ® 
it is intuitively clear that we will only see an ef- 
fect of the dynamical quarks on the microscopic 
spectral correlations if the quark mass is also of 
this size. For convenience, we rescale the quark 
mass by the same factor as the eigenvalues and 
define [i = mVY,. For /i > 1, we should simply 
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obtain results which are identical to the quenched 
approximation, since in this limit the fermion de- 
terminants do not affect the small eigenvalues. To 
observe results which are different from those ob- 
tained in the quenched approximation, we require 

In the following, we summarize recent results 
obtained in Ref. @. Analytical RMT results for 
the microscopic spectral density in the presence 
of massive dynamical quarks are currently only 
known for the unitary (UE) and chiral unitary 
(chUE) ensemble jF3|. For the chSE, results are 
known in the chiral limit. We have p^,p|,pT| 

p s (z) = z[Jl(2z)- J Q+ i(2z)J Q _ 1 (2z)] 

1 f 2z 
--J a (2z) / dtJ a (t) , (7) 
z Jo 

where a = Nf + 2\v\ and J denotes the Bessel 
function. If a is zero or an odd integer, RMT 
results for the distribution of the smallest eigen- 
value are known analytically fl5| , |i"6|| . For even 
a ^ (this is relevant in our case), they can 
be obtained as an expansion in zonal polyno- 
mials @|0i- For fi ^ 0, the RMT results 
for p s (z) and -P(A m i n ) can be obtained numeri- 
cally by constructing skew-orthogonal polynomi- 
als which obey orthogonality relations determined 
by a weight function involving the fermion deter- 
minants [ pTp^ ]. To avoid cancellation problems, 
a multi-precision package was used |L8| . 

Some of our numerical results are shown in 
Figs. H and [I| In this first exploratory study, 
we chose values of (5 in the strong-coupling region 
since there one does not need very large lattices to 
obtain agreement with RMT, see the end of Sec. || 
and Sec. [| The RMT results were computed with 
v = for the reason discussed in Sec. f|. Again, 
there is only one parameter, VS, which sets the 
energy scale and is determined by the data from 
the Banks-Casher formula. We have also plot- 
ted the RMT results for p s (z) and -P(A m ; n ) in the 
quenched approximation and in the chiral limit, 
respectively. The data should agree with these 
two curves in the limits fi ^> 1 and \x = 0, re- 
spectively. For [i ~ 0(1), the data should lie 
somewhere in between these two limiting curves. 
This is indeed observed in the figures. 
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Figure 3. Distribution of the smallest eigenvalue 
(top) and microscopic spectral density (bottom) 
for the simulation parameters indicated above 
the figures. The histograms represent the lat- 
tice data. The solid curves are the RMT re- 
sults computed with the appropriate value of the 
rescaled quark mass p. The dashed and the dot- 
ted curves are the RMT results computed in the 
quenched approximation and in the chiral limit, 
respectively. 



An important remark is in order. We are com- 
paring lattice data computed with n copies of 
staggered flavors (corresponding to Nf — 4n fla- 
vors in the continuum limit) with RMT results 
computed for Nf flavors, where Nf — Nf/2. This 
is not an ad-hoc prescription but justified by the 
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Figure 4. Same as Fig. || but for different simula- 
tion parameters (indicated above the figures). 



combination of the following two reasons. First, 
the U(4)xU(4) symmetry of the continuum ac- 
tion (for N f = 4) is broken to U(l)xU(l) at 
finite lattice spacing a. This would suggest to 
use Nf = Nf/4 in the RMT results. Second, all 
eigenvalues of D arc doubly degenerate in SU(2). 
However, in Eq. (|j) the eigenvalues are assumed 
to be non-degenerate ||. This would suggest to 
use Nf — 2Nf in the RMT results. Combining 
the two factors, we conclude that we have to use 
Nf = \ ■ 2 ■ N f = Nf/2 in the RMT results. 

Note that in the regime fi ~ 0(1), the quark 
mass is very small, and lattice simulations are 
very demanding. Testing the RMT predictions 



required a substantial numerical effort. However, 
the point is that, once the validity of the RMT 
description is established, the availability of an- 
alytical results for very small quark mass should 
facilitate extrapolations to the chiral limit which 
are otherwise difficult to perform on the lattice. 

6. Crossover to non-universal behavior 

As indicated above, the domain of validity of 
the random matrix description is finite. After all, 
QCD is a complicated theory and not just a ran- 
dom matrix model. If one wants to make prac- 
tical use of the RMT results, one should know 
quantitatively in which energy range they apply. 
In condensed matter physics, the energy up to 
which the RMT description is valid is called the 
Thouless energy, E c . Essentially, E c is the in- 
verse of the diffusion time of an electron through 
a disordered mesoscopic sample. It behaves like 
E c ~ L~ 2 , where L is the length of the sam- 
ple. Recently, some progress was made in iden- 
tifying the equivalent of the Thouless energy in 
QCD [^9|-p3[. Again, we shall concentrate on the 
low- lying eigenvalues. 

The starting point is the condition for the ap- 
plicability of the Leutwyler-Smilga effective the- 
ory, 



L<-L, 



(8) 



which also sets the domain of validity of the RMT 
description. From the Gell-Mann-Oakes-Renner 
relation, we have 



mlfl = 2m£ 



(9) 



Here, m is to be regarded as a valence quark mass 
which sets the energy scale below which RMT ap- 
plies. Denoting the QCD equivalent of the Thou- 
less energy by Armt, one obtains l20|j2l[| 



A 



RMT 



f 2 

SL 2 



(10) 



To have a dimensionless quantity, it is useful to 
divide Armt by the mean level spacing at A = 0, 
A = l/p(0) = 7r/(V£). This yields 



Armt/ A ~ —f%L 2 

TT 



(11) 
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This is the prediction which we want to test. It 
was already checked in the instanton liquid model 
and confirmed qualitatively [|T) . Since we have a 
large number of complete spectra at our disposal, 
we are in a position to perform a comprehensive 
quantitative test against lattice data. For simplic- 
ity, we only use the quenched data. In the follow- 
ing, we summarize results obtained in Ref. 22 1. 

A convenient quantity to consider in this re- 
spect is the disconnected chiral susceptibility. It 
is defined in terms of the Dirac eigenvalues by 



disc 
A lattice 



N 



1 



\kl~ll ^ k + m )(*^ + m ) 
, 2 

1 1 \ 




f-J iX k + m 

k—1 i 



(12) 



where N — Va~ 4 is the number of lattice sites 
and m is the valence quark mass. Going over to 
the microscopic limit, i.e., rescaling all energies 
by VTi and defining u — mVT,, the susceptibil- 
ity can be expressed in terms of the microscopic 
spectral one- and two-point functions. Inserting 
the RMT results of the chSE for these quantities 
and performing some tedious algebra, we obtain 
the RMT prediction 



Xrmt = 4u 2 



l r i 

2 ] 



dss z K Q (2su) / dtl (2stu) 
o Jo 

|s(l - t 2 ) - 8stI (2stu)K (2su) 
+AK Q (2u) [I (2su) + tl (2stu)]\ 



-Au 2 Kl(2u) 



ds Iq(2su) 



(13) 



where / and K denote modified Bessel functions. 
Note that in going from Eq. (P) to Eq. (|l3|), x disc 
has been rescaled by 1/ ( NYr) to eliminate the 
dependence on (3. The RMT result should agree 
with the lattice result up to some limiting value 
wrmt- To identify this value, it is useful to define 
the ratio 



ratio 



v disc 
Alatticc 



^ disc 
ARMT 



...disc 
ARMT 



(14) 



Deviations of this ratio from zero indicate the 
breakdown of the RMT description. The ratio 
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Figure 5. The ratio defined in Eq. ( |l4| ) versus 
u for the data obtained on a 10 4 lattice using 
(3 = 2.0. 



is plotted in Fig. || for one particular data set. 
The deviations of the ratio from zero for very 
small values of u are artifacts of finite statistics 
which are explained in Ref. |p2| . We are inter- 
ested in the deviations from zero which set in 
at u m 5 ~ 10. Let us denote this value of 
u by mrmt- It is related to Armt by urmt = 
XrmtVT, = 7tArmt/A. 

According to Eq. ( JTl] ) , wrmt should scale with 
L 2 . This prediction is tested in Fig. || where we 
have plotted data for four different lattice sizes 
at constant /3 versus u/L 2 . All data fall on the 
same curve, except for the deviations at small 
u which are due to the finiteness of our statisti- 
cal ensembles. This confirms the prediction that 
wrmt scales with L 2 . 

It remains to test the predicted scaling with 
f 2 . In order to determine /„•, we use a result of 
Ref. |24f| where it was found that f 2 w 3.4S in lat- 
tice units in the range of (3 we consider. This sug- 
gests to plot the ratio of Eq. (|l4| ) versus w/(£L 2 ) 
which we have done in Fig. [7j Again, the data 
fall approximately on the same curve, in particu- 
lar the data for (3 = 2.0 and (3 = 2.4. (One can 
only expect an agreement on the level to which 
the relation between f 2 and S is known.) Putting 
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Figure 6. The ratio of Eq. ( [l4| ) versus u/L 2 for 
constant (3 — 2.0 and four different lattice sizes 
V = L A . 
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A L= 8 , jj= 2.20 
■ L=6,/3=2.20 




u/(EL 2 ) 



Figure 7. The ratio of Eq. (|TJ) versus u/(SL 2 ) 
for four different parameter sets. 



in the numbers, we obtain 

A R mt/A « 0.3/ 2 L 2 (15) 

which is in good agreement with the prediction of 
Eq. (H|). 

Another interesting quantity is the number 
variance in an interval / = [0,5], defined by 
£ 2 (/) = <(jV(J) - (N(I))) 2 ). (This £ 2 should 
not be confused with the absolute value of the 
chiral condensate.) Here, N(I) is the number of 
eigenvalues in J, and (• • •) denotes an ensemble 
average. The number variance can be computed 
in RMT, and the lattice data agree with the RMT 
result up to some limiting value of S which we 
will denote by Srmt- Considering a number of 
different lattice sizes and /3-values we found that 
Srmt scales with f 2 and L 2 as expected. Quan- 
titatively, 

Srmt « (0-3 ~ 0.7)f 2 L 2 (16) 

which is consistent with urmt- 

We conclude this section with the statement 
that the domain of validity of the RMT descrip- 
tion of the low-lying Dirac eigenvalues in QCD 
is now known quantitatively. The implications of 
this result will be discussed in the next section. 



7. Summary and outlook 

We now have solid numerical evidence that the 
distribution of the low-lying Dirac eigenvalues is 
universal and described by RMT results, both in 
the quenched approximation and with dynamical 
fermions. Furthermore, we have a quantitative 
criterion for the domain of validity of the random 
matrix description which was confirmed by lattice 
simulations. How can we make practical use of 
this knowledge? 

One point concerns extrapolations to limits 
that are difficult to take on the lattice. The 
first example is the thermodynamic limit. Since 
we are dealing with finite-volume effects, we can 
make analytical statements about how this limit 
is approached. This was already demonstrated in 
Ref. ||. The second example is the chiral limit. 
We have RMT results for the microscopic spectral 
quantities in the presence of dynamical quarks 
with arbitrarily small masses. For quantities that 
are sensitive to the small eigenvalues, these re- 
sults should provide guidance for extrapolations 
to the chiral limit. Third, we have the continuum 
limit. It is less clear in what way the random 
matrix approach might help in approaching this 
limit. One observation is that as a — ► 0, there 
should be a transition in the effective number of 
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flavors used in the RMT results with dynamical 
fermions, from Nf/2 to 2Nf. Presumably, a has 
to be very small to see such a transition. 

As far as topology is concerned, RMT can de- 
termine the microscopic spectral quantities only 
for fixed topological charge v. The overall result 
for some quantity will be a weighted average over 
all topological sectors, and RMT cannot predict 
the weights. So far, all data were consistent with 
v = in the RMT results since the latter are 
only sensitive to exact zero modes of the Dirac 
operator. It would be very interesting to identify 
the would-be zero modes. For this, one needs the 
eigenvectors corresponding to the small eigenval- 
ues. We are currently trying to compute these as 
well. 

Last but not least, we feel that the analyti- 
cal knowledge of the distribution of the smallest 
eigenvalues may have an impact on fermion algo- 
rithms. This is suggested by the fact that the 
magnitude of the small eigenvalues determines 
the performance of the fermion algorithm. While 
RMT cannot say anything about individual con- 
figurations, it predicts the statistical distribution 
of the small eigenvalues. This knowledge may be 
useful to construct more effective algorithms. 
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